Phys. Rev. D 82, 083006 (2010) 



arXi v:1001. 1939 



Effective cosmological constant from TeV— scale physics 



F.R. Klinkhame 

Institute for Theoretical Physics, University of Karlsruhe, 
Karlsruhe Institute of Technology, 
CN ' 76128 Karlsruhe, Germany 



Abstract 



o 
O 
o . 

It has been suggested previously that the observed cosmological constant A corresponds to the 
remnant vacuum energy density of dynamical processes taking place at a cosmic age set by the mass 
scale M ~ E cw of ultramassive particles with electroweak interactions. Here, a simple modification 
of the nondissipative dynamic equations of g-theory is presented, which produces a remnant vacuum 
energy density (effective cosmological constant) of the correct order of magnitude. Combined with 
the observed value of A, a first estimate of the required value of the energy scale E ew ranges from 3 
to 9 TeV, depending on the number of species of ultramassive particles and assuming a dissipative 
Q\ ' coupling constant of order unity. If correct, this estimate implies the existence of new TeV-scale 
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I. INTRODUCTION 



It has been argued by Arkani-Hamed et al. [l( that two fundamental energy scales, the 
electroweak scale E ew ~ 1 TeV and the gravitational scale -Epianck ~ 10 15 TeV, suffice to 
explain the triple cosmic coincidence puzzle: why are the orders of magnitude of the energy 
densities of vacuum, matter, and radiation approximately the same in the present Universe? 
For this explanation to work, the parametric form of the effective cosmological constant 
(remnant vacuum energy density) must be 

A = PV,remnant ~ {(E cw ) 2 /^Planck)' ~ (1(T 3 eV) 4 . (1.1) 

If true, formula ( II -ip would be a remarkable explanation of the measured value from ob- 
servational cosmology, which appears to be of order 10~ 29 gem -3 ~ KT 11 eV 4 (setting 
h — c — 1 and referring to, e.g., Refs. 2|-|4j and other references therein). 1 However, (II. ip 
was not derived convincingly in Ref. [1], as an unknown adjustment mechanism needed to 
be invoked. 

Subsequently, Volovik and the present author realized jsl that, in the framework of 
g-theory, there is the possibility of generating a vacuum energy density precisely of the 
form (II. ip . Here, g-theory is a particular approach j3] to solving the first cosmological 
constant problem (CCP1): why is |A| <C (-Epianck) 4 ? The original references on the statics 
and dynamics of g-theory are |8| and [9j, respectively. The second cosmological constant 
problem (CCP2) is the question addressed here, namely, the actual order of magnitude of 
A, if indeed nonzero. 

The positive remnant vacuum energy density obtained in Ref. relied erucially on 
Eq. (4.1) of that article. That particular equation was taken to describe the quantum- 
dissipative effects of the vacuum energy density, but was, in the end, purely hypothetical 
and disconnected from the previous g-theory discussion. 

The question arises if it is at all possible to modify the previous g-theory equations 
in such a way as to effectively recover the results of Ref. [6j. The answer is that this is 
indeed possible, even though the required modifications are quite subtle. The scope of the 
present article is restricted to finding these appropriate phenomenological equations, rather 
than establishing the relevant microscopic processes of the underlying theory. 

In a way, it can be said that this whole article is about the proportionality constant 
implicit in ( II. ip and that, within the framework of g-theory, the article gives an existence 
proof for a set of dynamical equations which produces a proportionality constant of order 
unity. 

As the previous discussion makes clear, the present article is a direct follow-up of Ref. j^], 
to which the reader is referred for the original motivation and detailed analysis. This article 



There have, of course, been many other explanations of the smallness of A by an appropriate ratio of 
energy scales (see, e.g., Sec. X of Ref. [5|]), but the relation (jl.ll) is special as it carries the ingredients to 
naturally give the correct orders of magnitude for the present matter and radiation energy densities 
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is, by necessity, rather technical and it may be helpful to give the reader a road map. The 
material of this article is organized as follows: 

Track 1: Sees. HQ ITTTAl UTTBl and [VI 

Track 2: Sees. IriTH UTTDl andHVJ 

Track 3: Sec. HTHland Appendix \M 

The basic idea and main results are presented in Track 1, the dimensionless differential 
equations and their numerical solution in Track 2, and a more detailed discussion and further 
refinement in Track 3 (the most realistic calculations are shown in the very last two figures 
and the very last table of Appendix [A}. In a first reading, it is possible to follow Track 1 
and to add the other Tracks later. 



II. THEORETICAL FRAMEWORK 

This section reviews the main ingredients of the type of theory considered in this article 
(see Refs. jg[ lof for details). The particular g-theory realization used involves the so-called 



4- form field strength [10|, [III]. Very briefly, the theory is defined over a four- dimensional 
Lorentzian spacetime manifold and employs a 4-form field strength F derived from a 3-form 
gauge field A. The corresponding rank-four tensor can always be written as 



F a ^s(x) = q(x) \/-g(x) 6 a /3 lS (x) , (2.1) 

with the Levi-Civita tensor density e a/ 3^s(x), the determinant of the metric g(x) = 
det(7 Q/ g(x), and the scalar field q(x). The crucial point is that this scalar field q(x) is 
nonfundamental, being built from the metric field g a /3(x) and the 3-form gauge field A(x), 
as will become clear shortly. This 3-form gauge field A(x) is considered to be one of the 
fields which characterize the quantum vacuum at the fundamental microscopic level (sj. 

The macroscopic effective action of the relevant high-energy fields (here, A) and the 
low-energy fields (here, g and ip) is taken to be of the following form 0|: 

S cS [A,g,ij}= [ d'x^(K(q)R[g} + e(q) + C M [ij,g}), (2.2a) 
F a p~t8 = V[aAp^s\ , (2.2b) 
q 2 ee -If^F*™, (2.2c) 

where the effective gravitational coupling parameter K is allowed to depend on q, R[g] is 
the Ricci curvature scalar obtained from the metric g a p, V Q denotes the standard covariant 
derivative, and the square bracket around spacetime indices stands for complete antisym- 
metrization. The energy density e(q) is assumed to be a generic function of q, that is, 
a function different from the simple quadratic \ q 2 corresponding to a Maxwell-type the- 
ory 10, 111) . The field ip in (12.2aj) stands for a generic low-energy matter field with a scalar 
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Lagrange density C M [tp,g], which, for simplicity, is assumed to be without explicit g-field 
dependence (the dependence on the metric arises from the covariant derivatives). 



Remark that the effective action (12.2aj) corresponds to a Brans-Dicke-type action [12] , but 
without kinetic term for the (nonfundamental) scalar q. For spacetime-independent q [that 
is, q(x) = q = const], the effective action (12.2aj) corresponds to the one of standard general 
relativity with a cosmological constant A = e(q) + A M , where A M refers to contributions to 
A from the matter Lagrange density C M . 

By taking variations of A a ^(x) and g a p{x) in the effective action (12.2a[) . generalized 
Maxwell and Einstein equations can be derived. The generalized Maxwell equation can 
be solved explicitly and the solution depends on a constant of integration fi. With the 
solution of the generalized Maxwell equation, the generalized Einstein equation reduces to 
the following field equation: 

2K (R a0 - g af3 R/2) = -2 (V a V^ - g afi □) K(q) + p v (q) g afi - , (2.3) 
where the combination 

Pv(q) = e{q)-nq (2.4) 

plays the role of the gravitating vacuum energy density rather than the single term e(q) 
appearing in the effective action (I2.2al) . Furthermore, there is an equation remaining from 
the particular solution of the generalized Maxwell equation, which reads 

d -^ + R d -f = 0. (2.5) 
dq dq 

The final equations f l2.3p -( l2~5l) can be specialized to the case of a spatially flat Friedmann- 
Robertson-Walker universe. The resulting cosmological equations have been studied in 
Refs. 0, 9] and it is the aim of the present article to find a modification of them which 
allows for the generation of a nonvanishing remnant vacuum energy density. 



III. VACUUM DYNAMICS IN A FLAT FRW UNIVERSE 
A. Basic idea 

Following Ref. [lj], assume the existence of ultramassive unstable particles (here, called 
'type la') with masses M of order E ew ~ 1 TeV and electroweak interactions. Consider a 
spatially flat Friedmann-Robertson-Walker (FRW) universe and assume the type-la par- 
ticles to be effectively in thermal equilibrium at early enough times. Then, the masses of 
these particles start to affect the Hubble expansion rate H{t) when the temperature drops 
to T ~ E ew , corresponding to a cosmic age t of order 

tew = -E'pianck/(-Eew) 2 , (3.1) 

in terms of the reduced Planck energy, 

^Planck = y/l/(S%G N ) w 2.44 x 10 18 GeV. (3.2) 
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Note that the definition (13. ip is motivated by the standard Friedmann equation 
H 2 ~ p/(-Epianck) 2 with H ~ 1/tew and p ~ T 4 ~ (E ew ) 4 . In the following, the Friedmann 
equation will be modified, but the order of magnitude ( 13. ip remains relevant. Throughout 
this article, natural units are used with H — c—1. 

Compared to the case of having only ultrarelativistic particles (these lighter particles 
are called 'type lb' and can be thought to have masses of order M/10), the change of 
the expansion rate from the ultramassive type-la particles can be modeled by a nonzero 
function KMi{t/t ew ). In fact, this function kmi can be written in terms of the standard 
equation-of-state (EOS) parameter wmi = PmxI ' Pmi as follows: 

«mi = 1 -3idmi , (3.3) 

which vanishes for ultrarelativistic particles (wmi = 1/3) and equals unity for pressure- 
less nonrelativistic particles (wmi = 0). By taking different unstable type-1 particles 
(a, b, c, . . .) it is possible to obtain an effective EOS function KMi(t/t ew ) which peaks at 
t = t ew . Here, however, a particular form of Kmi will simply be assumed. 

The main conditions on this assumed EOS function kmi are that it peaks at t — t ew and 
is nonzero only in a finite range around the maximum (this last condition is not essential 
but simplifies the discussion). Specifically, the conditions are taken to be: 

KMi(t/tew) < «Afi(l) for t^t ew , (3.4a) 



with < t s tart < t ew and t ew < t en d < oo, having set t = for the big bang where H(t) 
diverges. The physical picture corresponding to (13 .4p is that ultramassive type-la particles 
are dominant at t ~ t ew , but, then, decay into lighter type-lb particles which are still 
ultrarelativistic for t not very much larger than t ew . The main goal is to study the effects of 
this prescribed EOS function Kmi, relegating the discussion of a more realistic EOS function 
to App. \M 

A flat FRW universe containing only type-1 particles with a prescribed EOS parameter 
(l3.3p - (l3.4p has a standard radiation-dominated Hubble expansion rate H(t) = (1/2) t' 1 for 
t < tstart and t > t en d. The expansion rate is changed, H(t) ^ (1/2) i -1 , for times t between 
t s tart and tend- The question, now, is what happens if this FRW universe also has a dynamical 
vacuum-energy-density component. 

For the theory outlined in Sec. [TTJ two results were obtained in Ref. pj. First, it was 
shown that there is an exact solution having pv(t) = in the radiation-dominated phase 
with KM\{t) = 0. Second, it was shown that the changed Hubble expansion from KMi(t) ^ 



kicks pv(t) away from zero. Specifically, the following behavior was established [6| 



which vanishes asymptotically as kmi drops to zero and the standard radiation-dominated 
expansion of the model universe resumes. At the moment of the kick, t ~ t ew , the vacuum 



K M l(t/tew) 7^ 
KMl(t/tew) = 



for t G (t s tart, tend) ? 
for t ^ (tstart, tend) , 



(3.4b) 

(3.4c) 




(3.5) 
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energy density (jH3J) is of order (t ew )~ 4 ~ ((E cw ) 2 / -Epianck) 4 ~ (E ew /E Planck ) A (E ew ) 4 , which 
is negligible compared to the matter energy density pu\ ~ (-^cw) 4 - The vacuum energy 
density pv(t), therefore, just responds to (is being kicked by) the Hubble expansion and 
does not affect the expansion substantially. 

The result (I3.5P has been obtained from the simplest time-reversible (nondissipative) 
version of g-theory, with field equations given by (12.31) — ( T2~5T) . It has been argued that 



quantum-dissipative effects (e.g., because of particle production in an expanding universe 13 



14|) may result in a freezing of the previous result (13. 5 j) to a constant nonzero value. 

As explained in Sec. HI the aim of this article is to find a suitable modification of the 
"classical" g-theory equations, which produces a finite remnant vacuum energy density. In 
the approach followed here, there are three changes: 

1. The matter energy-conservation equation is modified to include appropriate particle- 
production effects operating at a cosmic age t ~ £ ew . 

2. The reduced Maxwell equation is modified, so as to match the standard Einstein 
equation of an FRW universe with a nonzero effective cosmological constant at later 
times. 

3. Different particle species are considered with ultramassive type-la particles first decay- 
ing into lighter type-lb particles, which, in turn, decay into massless type-2 particles. 

The first two modifications are essential for the generation of a nonvanishing remnant vacuum 
energy density. (The first modification has already been discussed in general terms in Sec. IV 
of Ref. jsl). The third modification allows for a possibly more realistic scenario, with type-1 
particles corresponding to new TeV-scale physics and type-2 particles corresponding to the 
standard model of elementary particle physics (see also App. lAl . 

The particle-production effects of point 1 above will be controlled by an effective coupling 
constant ( > and a particular type of dissipation function 7(£/t ew ; () > 0. The reason for 
calling 7 a "dissipation" function will become clear in Sec. IIIIEI The main conditions on 
this function are as follows: 

7(0; C) = 1, (3.6a) 
Vt>tf r eeze: l(t/t cw ; () = 0, (3.6b) 

lim (tfreeze)" 1 = , (3.6c) 

with a particular time tf reeze that is of order t ew for ( ~ 1 and approaches infinity for ( — > 0. 

The coupling constant ( and the corresponding function 7 are purely phenomenological. 
The vanishing of 7 for large enough times will be seen to have two effects: first, to freeze the 
"classical" value (13.51) and, second, to switch to a standard FRW expansion with relativistic 
matter and a tiny value (11. ip for the remnant vacuum energy density. Further discussion 
of ( and 7 will appear in the next subsection, after the differential equations have been 
presented. 
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B. Modified ODEs 



For a spatially flat FRW universe and the 4-form realization of g-theory with a variable 
gravitational coupling parameter K(q), the cosmological differential equations have been 
derived in Ref. [9| and were already mentioned in Sec. [TTJ The basic idea of the proposed 
modification of these ordinary differential equations (ODEs) has been discussed in the pre- 
vious subsection. Specifically, the modified ODEs are given by: 



dp 



Ml 



dt 



dK ( dH 9 \ r . . ,-, dpv r , i ,i 1 dK 

6 W U + 2H J = Ht/t - )] l{t/t - )] K~q 2p ^ ^ 

+ [4 - K M1 (t/t ew )] Hp M1 = <4 (*%L) - ^ [1 - j(t/t ew )] p M1 , (3.7b) 

7(t/t cw ) dt V dq J t ew 

^Ml+4Hp M2 = +^ [1 - 7 (*/*ew)] Pmi , (3.7c) 

Lib 6g w 

6 (h^^ + KH 2 ) = p v + p M1 + p M2 , (3.7d) 



where only the arguments of the functions Kmi and 7 have been shown explicitly. The four 
equations in ( 13. 7ft . going from the top to the bottom, can be recognized as modified versions 
of the reduced Maxwell equation (12.51) . the two matter energy- conservation equations, and 
the Friedmann equation [the standard Friedmann equation is recovered for K = K(q ) = 
1/(16ttG n ) = const]. For ( = and 7(f) = 1, the ODEs fl32p correspond to Eqs. (4.12abc) 
of Ref. [91], supplemented by an equation for the adiabatic evolution of pu2- 

The modified ODEs for ( > have a dependence on an external time scale, here taken to 
be t ew from (13.11) . The t ew dependence enters implicitly through the EOS function k m1 > 
and the dissipation function 7 > discussed in the previous subsection and explicitly through 
the matter energy-exchange terms proportional to Xu/tew Note that the dimension of the 
g-field in Eqs. (13. 7p is irrelevant, which concords with the fact that this g-field may be 
realized in different ways jsj. 

The particular modification ( 13.71) of the cosmological ODEs from classical g-theory has 
two main ingredients: first, the function , y(t/t ew ) with characteristics H 3 . 6 1) and, second, the 
presence of a finite coupling constant (, 

C = 0(1), (3.8) 

which enters directly on the right-hand side of (I3.7bj) and indirectly via condition (I3.6cl) for 
the dissipation function 7. A finite value for the remnant vacuum energy density from the 
dynamical ODEs (13 .7p requires both ( > and ( < 00, as will be explained in Sec. IIII El 

The rather simple structure of (13.71) . combined with conditions (13. 4p and (13. 6p . will be 
seen to allow for the generation of a nonzero remnant vacuum energy density. 2 A detailed 



2 A somewhat more general modification of the ODEs has 2pv in the last term on the right-hand side of 
(|3.7a[) replaced by [2py + KMi{t/t cw ) pmi/2], but, for the case considered in this section and SecHV] the 
results are essentially unchanged, because KMi(t/tcw) — for t > t cn d > t cw according to (|3.4[) . 
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discussion of the modified ODEs is postponed until Sec. IHI El after these equations have 
been established in dimensionless form. 

Before embarking on this technical enterprise, it may be useful to recapitulate the basic 
assumptions. The first assumption is the existence of a particular type of vacuum variable 
q(x), namely, a variable which corresponds to a conserved relativistic quantity go i n flat 
Minkowski spacetime. Such a variable q(x) provides a possible solution of the main cos- 
mological constant problem (CCP1) by explaining why A/£p lanck is naturally zero in the 
equilibrium state. 3 This vacuum variable q(x) is taken to have an effective action of the form 
of (12.21) . where, in particular, the gravitational coupling constant K may carry a dependence 
on q. 

The second assumption is that the field equations from the effective action (12. 2p . special- 
ized to a spatially flat FRW universe, are modified by the introduction of terms involving the 
coupling constant (. The crucial term is the first one on the right-hand side of (]3.7b|) , whose 
physical motivation is that it reproduces the dissipative behavior suggested in Ref. jg] (this 
behavior is analogous to that of bulk viscosity in compressible material fluids 15|]). The 
coupling constant ( and the corresponding function 7 are purely phenomenological. As 
mentioned in Sec. [U ultimately ( and 7 (or appropriate generalizations) need to be derived 
from the underlying microscopic theory, but that task lies outside the scope of the present 
article. 

Clearly, the first assumption is better motivated than the second. But the second assump- 
tion may (or may not) gain in credibility depending on the success (or not) of producing a 
reasonable remnant vacuum energy density and predicting new TeV-scale physics. 



C. Ansatze and dimensionless variables 

Following Refs. 0, ^ , take quadratic and linear Ansatze for the vacuum energy density 
and the gravitational coupling parameter: 

1 2 

Pv(q) = 2 (l ~ <?o) , (3.9a) 

K(q) = X -q. (3.9b) 
These Ansatze imply the following equilibrium value for the g-field of theory (12.21) : 

q = l/(8nG N ) = (Planck) 2 , (3.10) 



3 The g-theory approach to CCP1 provides only a possible solution, because it is not known for sure that 
the underlying microscopic theory does contain an appropriate g-type field. In addition, there remain 
other equally fundamental (perhaps related) questions, such as the nature of gravity and the origin of 
spacetime. The goal of the present article is relatively modest: to explore, in the framework of q-theory, a 
possible connection between the observed value of the effective cosmological constant and new TeV-scale 
physics. 



s 



where -Epianck is the energy scale from (13.21) . In this article, q is considered to be realized 
by a 4- form field strength F with mass dimension 2. With a proportionality constant in 
( 13.9bj) of order unity, the natural scale of the vacuum variable q is then of order (-Epianck) 2 - 



However, the energy scale of q in the cosmological ODEs ( 13. 7p is, in principle, arbitrary, as 
noted already in the second paragraph of Sec. IIII Bl 

Next, recall the time scale t ew defined by (13.11) . which corresponds to the age of the 
Universe at a temperature of order E ew . For later use, also define the following number 
characterizing the hierarchy of energy densities: 

£ = (-Eplanck/^ew) 4 , (3-H) 

which is of order 10 60 for E ew ~ 1 TeV. For such a large value of £, the cosmic time t ew 
considered in this article is large compared to the Planck time, t ew = a/£ (£pianck) _1 - In 
addition, the relation (t ew ) 2 Qo = £ can be seen to hold, which will be used later for the 
derivation of the dimensionless ODEs. 

With £ ew and £, the following dimensionless variables can be defined for the cosmic time, 
the Hubble expansion rate, the energy densities, and the q shift away from equilibrium: 

r = (tew)" 1 *, h = t cw H, (3.12a) 

r v = (tew) A pv, r Mn = £ -1 (*ew) 4 PMn , (3.12b) 
x = £(?/<Zo-l) =£y, (3.12c) 

where n stands for the matter-species label (n = 1, 2) and y is the variable used previously 
in Refs. Observe that pu n has been rescaled by an extra factor l/£ but py not. 

At this moment, it is appropriate to give explicit examples for the EOS function kmx and 
the dissipation function 7 discussed in Sec. IIII Al With central value r c and total width Ar, 
define the auxiliary variable o = 2 (r — r c )/Ar and take the EOS function to be given by 

/ k c sin 2 [(tt/2) (1 + a 2 )] for - 1 < <r < 1 , 
{ otherwise . 

In the main part of this article, k c is set to 1 and the dynamic strength of the kick is 
controlled by the value of the initial energy- density ratio Pm\I Pmi ( see App. [A] for a more 
realistic EOS function). 

Turning to the dissipation function 7, introduce the basic time scale Too, define 

Tfreeze = (1 + 1/0 Too, (3.14a) 

and take the function to be given by 

( r ) _ / C0S2 K 71 / 2 ) T / r freeze] for < T < T beeze , 

1[T) ~ \ otherwise . ' 

The two functions used will also be shown in the plots of the numerical results later on. 
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D. Dimensionless ODEs 



Take, now, the Ansdtze (I3.9ap -f l3.9bj) and assume small deviations of q away from the 



equilibrium value q , i.e., |g/ g — 1| <C 1. Then, the ODEs (13.71) reduce to the following four 
equations for the four dimensionless variables h(r), tmi(t), Tm2{t), and x(t) from A3. 12[) : 

3(h + 2h 2 ) = 7 x + (1 - 7) £ -1 x 2 , (3.15a) 

Tmi + (4 - K M i) hr M1 = -(C/t) [&] ~ A12 (1 - 7) r Mi , (3.15b) 

f M2 + 4:hr M2 = +A12 (1 - 7) r M i , (3.15c) 

3 (r 1 h x + h 2 ) = r 1 x 2 /2 + r M1 + r M2 , (3.15d) 

where the overdot stands for differentiation with respect to r. Henceforth, the functions 
k>mi(t) and 7(7-) are considered to be given by the explicit expressions ( 13 . 13[) and (j3.14j) . 
For the numerical calculation, the factor [x] on the right-hand side of f!3.15bj) is to be replaced 



by the appropriate expression for x obtained from (I3.15dj) . Recall, furthermore, that the 



dimensionless vacuum energy density ry is given by x 2 /2 and that, according to (13.12bj) . the 
dimensionless matter energy densities Tmi and tm2 include an extra numerical factor 
compared to r v . 

For later use, also the ODEs for the special case = are needed. From (I3.15p . the 
following system of equations can be derived for = 0, with three ODEs: 

3 (h + 2h 2 ) = 71, (3.16a) 

x = -C 1 hj (2jx - Kmi [3/i 2 - r M2 ] ) , (3.16b) 

r M2 + 4/ir M2 = +A12 (1 - 7) [3/i 2 -r M2 ] , (3.16c) 
and a single algebraic equation: 

3h 2 = r M \ + r M2 . (3.17) 

The derivation of (13 . 16b|) proceeds in three steps: first, take the time derivative of (j3.17|) : 
second, use (I3.16aj) to eliminate h in the resulting expression for {tmi + ^a/2); third, use 
the final expression for (r M1 + f M2 ) in the sum of the two Eqs. f!3.15bj) and (13. 15cj) to get 



( I3.16bj) . All in all, the ^ 1 = equations consist of three ODEs f l3.16ap - fl3.16cp for three 



variables h(r), x{r), and ^(t), with the energy density tmi following from the Friedmann 
Eq. (I3.17p . These ODEs will be used in Sec. HVBI to get the value of the remnant vacuum 
energy density ry for £ — > 00. 

Purely mathematically, there is another special case to consider for the ODEs (I3.15P as 
they stand, namely, the case £ = 0. From f)3.15ap and (I3.15d|) . together with the proper 



boundary condition x(0) = 0, the £ = solution can be seen to have x(t) = 0, which implies 
ry(r) = 0. But, as said, this solution is not directly relevant for the physical situation 
considered. 
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E. Additional remarks 



A few remarks may be helpful to better understand the proposed ODEs, given by (13. 7p 
in the general form or by (13.151) in the specialized and dimensionless form. 

First, note that the previous results [6| on the dynamics at r ~ 1 are readily recovered. 
From (I3.15P for ( = 0, j(t) = 1, and £ Tmi 3> £ ^M2 S> ry, one immediately obtains at r ~ 1 

x(r) ~ - K M1 (r)h(r) 2 , (3.18a) 
r V (r) ~ I K M1 {r f h{r)\ (3.18b) 



which corresponds to Eqs. (3.1a) and (3.4a) of Ref. [6] apart from a trivial rescaling. With 
^mi(t) 2 M r ) 4 dropping to zero rapidly for large enough r, there is no sizable remnant vacuum 
energy density, at least, according to the unmodified ODEs given by (I3.15P for ( = and 
7 = 1. 

For nonzero (, however, the same approximations in the ODEs (13.151) give the following 
dissipation- type equation at r ~ 1: 

x(t) ~ - 7 diss(r) [j{t)x{t) - (3/2)K M1 (r)h(r) 2 ] , (3.19a) 

7 diss(r) = 2C 1 /i(r) 7 (r), (3.19b) 

whose derivation parallels the one of ( I3.16bl) in the previous subsection. Equation ( 13. 19aj) 
with boundary condition x(0) = is the analogue of the crucial relation (4.1) of Ref. [6| that 
allows for a positive remnant vacuum energy density as discussed in Sec. IV of that article. 
Section IV of Ref. contains, in fact, the analytic solution of (I3.19al) for given functions 
7(t), Kmi(t), and h{r). For completeness, the dimensionful quantity corresponding to 
7diss(V) is given by r diss (t) = 2 C^ 1 H(t) j(t/t ew ). 

The dissipative ODE (13.191) and its analytic solution [6| make clear that a finite remnant 
vacuum energy density with lim r _ >00 x(r) ^ requires both ( > and ( < oo. Indeed, for 
the case 7di SS — > oo (or ( — >• 0), the solution x{r) follows kmiM h(r) 2 which drops to zero 
rapidly for large times and, for the case 7di SS — > (or ( — >■ oo), the solution x(r) simply 
remains at the initial value, x(t) = x(0) = 0. 

Second, it appears essential that the ODEs ( 13. 71) and ( 13.15P are singular, with the coeffi- 
cients of the first terms on the right-hand sides of ( I3.7bl) and ( I3.15b[) diverging for r > Tf ree ze, 
because of the 7 condition (I3.6bl) . In fact, the divergence of the coefficient Cllij) f° r T above 
rf reeZ e forces (q/qo — 1) to be strictly constant. It could very well be that the exact vanish- 
ing of 7 for r > rf re e Z e in the cosmological context traces back to the existence of an energy 
threshold in the relevant particle reaction process. As noted in Ref. 0] , the energies involved 
are tiny (of the order of meV), so that only sufficiently light neutrinos and gravitons can be 
expected to play a role. The first term on the right-hand side of (I3.7bl) . as it stands, is not 



simply proportional to R 2 as for the well-known Zeldovich-Starobinsky result 131 ]. but does 
involve R via its time derivative, as follows by use of (I3.7aj) . 
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Third, the discussion of the two previous remarks suggests that the value of the remnant 
vacuum energy density can be at most of the order of the maximum possible "classical" value 
(i.e., the ry peak from the nondissipative ( = equations). The idea is that, in general, 
dissipation leads to reduction of the produced energy rather than enhancement. Specifically, 
the conjectured inequality is 



?V( r free z e) 



< 



l^Mi(r) 2 {h(r) 2 -r M2 (r)/3Y 



(3.20) 



which is based on the analytic result f!3.18bj) with h 4 on the right-hand side replaced by 



{h 2 — rM2/3) 2 , as suggested by ( 13.16bl) . It remains to sharpen the approximate upper bound 
(I3.20p and to determine the corresponding conditions. 

Fourth, having a constant nonzero value of x oc (q/qo — 1) does not automatically allow 
for a standard de-Sitter universe, as the original ODE f)3.15ap [with 7 = 1] and the ODE 
f!3.15dl) are inconsistent for x = h = 0, x < £, and Tmi = ?"M2 = 0. However, the modified 
Eq. (13.15aj) [with 7(7-) = for r > Tf reeze ] has been designed to match the corresponding 
Einstein equation of a standard flat FRW model with ultrarelativistic matter and constant 
vacuum energy density, which asymptotically approaches a de-Sitter universe. 4 This par- 
ticular modification also makes clear that there must be more than just energy exchange 
between the vacuum and matter sectors. Rather, there must be a type of modulated inter- 
action between the vacuum field and the nonstandard gravitational field, which can be seen 
as follows. Multiply (l3.7aj) by K dq/dK to get a modified FRW-Einstein equation, 

Gk(^ + 2H^j = [ T (t/tew)] T(V*ew)] 2p v , (3.21) 

where the nonstandard term K dpy/dK is switched off for large enough cosmic times by the 
factor j(t/t ew ) going to zero. 

Fifth, the ODEs (13.151) . for given functions kmi( t ) an d 7(7") and fixed coupling constants 
( and A12, contain one last free parameter, the hierarchy parameter £ defined by (13.111) . 
Heuristically, it is to be expected that the precise value of £ ^> 1 does not affect the resulting 
value ry (Tf r ee Z e)- But, as £ fixes the ratio pm/ pv at t ~ £ ew , it does affect the later (standard) 



4 The assumption, here, is that other contributions to the vacuum energy density generated at times later 
than t cw would be self-adjusted away by appropriate g-type fields 0]. The prime example is the quantum- 
chromodynamics (QCD) vacuum energy density of order (100 MeV) 4 , which is expected to appear during 
the cosmological QCD transition at T ~ 100 MeV. This huge contribution to the vacuum energy density 



pv has been shown to self-adjust to zero lsll2(j as long as there is no term proportional to |-ff|£Q CD 



contributing to py ■ If there is such a nonanalytic term (cf. the discussion in Ref. [18|), then the final value 
of A could be a combination of electroweak and QCD effects. However, the experimental signatures of the 
electroweak vacuum energy density (effectively, a ACDM model, as explained in the next footnote) and 
the QCD vacuum energy density (an f(R) modified-gravity model) are different, in particular as regards 



the effective EOS parameter wx discussed in Ref. 17|. But, for the moment, there is no definitive proof 



that the required nonanalytic term occurs in four-dimensional QCD. 
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evolution of the model universe and, in fact, determines |l| the cosmic time t acce i at which 
the matter energy density drops below that of the constant vacuum energy density, t accc i ~ 
\/£ t cw . More precisely, the onset of acceleration a/a > [having defined h = a/a in terms 
of the scale factor a(r)] occurs at the energy- density ratio pv / Putot = v, with v = 1 for 
relativistic matter (wm2 = 1/3 and km2 = 0) or v = 1/2 for nonrelativistic matter (wm2 = 
and km2 = !)• The model considered in the present article has v — 1, but can easily be 
adapted to give the value v = 1/2 which is more realistic. 5 



IV. NUMERICAL SOLUTION 



A. Numerical results for £ = 10 2 

The mathematical parameter £ entering the ODEs (I3.15P is, first, considered to have the 
moderately large value of 100. The boundary conditions are taken from the epoch before 
the electroweak kick, when there was a standard radiation-dominated flat FRW universe. 
With the onset of the electroweak kick given by r start from (13.4bl) and fl3.13p . the following 
boundary conditions on h, x, and fMtot = r Mi + r M2 hold at a time r = r min < r sta rt < t cw : 

h(r min ) = 1/2 (r min )- 1 , (4.1a) 

x(r min ) = 0, (4.1b) 

^Mtot(^min) = 3 [h(r min )] 2 . (4.1c) 

This leaves only the initial ratio [rjvn/rM2](7"min) = rMi(Tmin)/?"M2(T m in) to be determined, 
which is, for the moment, taken to be 1 (other initial ratios will be discussed shortly). 

The corresponding numerical solutions of the ODEs (I3.15P are shown in Figs. [THSJ 
Figure [1] illustrates the fact that the standard (nondissipative) dynamic equations for ( = 
and 7 = 1 do not produce a constant positive remnant vacuum energy density from the 
electroweak kick [the oscillatory effects in ry(r) are suppressed for larger values of the hi- 
erarchy parameter £; see the first figure called in Sec. UVBj . However, as shown by Fig. 121 
the modified (dissipative) dynamic equations for ( = 0(1) do produce a sizable remnant 
vacuum energy density. 

The subsequent evolution of Fig. [5] is shown in Fig. [31 The content of this model universe 
for t > 3 is given by a constant vacuum energy density (effective cosmological constant) and 
two species of matter, with massive type-1 particles playing a role for the generation of the 



5 The adapted model contains an additional EOS function KM2(t/t eq ), which is a smoothed step function 
running from to 1 as the cosmic time t increases and which has a half-way time t eq where km2(1) = 1/2. 
This matter-radiation-equality time t eq has a parametric form a 3 t ew , with a the electromagnetic fine- 
structure constant (see Ref. [1( for further details). The adapted model then has an accelerationphase for 
Pv/(Pmi + PM2) ~ Pv I Pmi > 1/2 and corresponds, for t ^> t cw , to a particular ACDM model [4|, 19]. As 
noted in Ref. [a], this EOS function Km2 can also be expected to perturb the vacuum energy density, but 
the magnitude involved is tiny compared to the one from the electroweak scale because iJ(i oq ) 4 <§C H(t ew ) 4 . 
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vacuum energy density during the electroweak epoch (r ~ 1) and ultimately decaying into 
massless type-2 particles. 

Similar results are obtained for initial ratios [rMi/ r M^[{ T tDm) ^ 1- Table Q] gives the func- 
tion values at cosmic time r = Tf reeze , in particular, the values for ry which, by construction, 
stay constant for later times r > Tf reeze . Remark that even for a relatively mild kick with 
initial ratio [r^i/r^] = 1/10 , the generated r v is still of order 10~ 3 . 

Returning to the boundary condition [r mi I f M2\ijmin) = 1, Figs. EH3] are seen to give 
a value 7v(Tf reezc ) w 0.04. By changing the model parameters and the model functions 
somewhat it is possible to get r v (r {reeze ) values in the range of 10 -3 to 1. But it appears 
impossible to get a remnant ry much larger than unity, which agrees with the conjectured 
upper bound (I3.20p . 



B. Numerical results for £ ^> 10 2 



The parameter £ has been defined in physical terms by (13.111) and its mathematical role 
for the solution of the ODEs (13. 15[) has already been discussed in the last paragraph of 
Sec. IIII El Here, numerical results are presented for large values of this parameter, ranging 
from £ = 10 4 to £ = oo. 

Numerical results for the standard nondissipative (( = 0) dynamic equations at £ = 10 4 
are given in Fig. HI which show reduced oscillatory effects of ry{r) compared to Fig. [I] and 
recover the smooth behavior of (13.18bj) . [Recall that the analytic approximation ( 13.18b|) 



was derived for a negligible type-2 matter energy density and a better approximation has 
/i 4 on the right-hand side replaced by {h 2 — rM2/3) 2 , as used already in (I3.20p .] Further 
numerical results for Q = 2 and £ = 10 4 confirm the expectation from Sec. IIII El that the 
generation of the remnant vacuum energy density at r ~ 1 is qualitatively the same as for 
£ = 10 2 (compare Fig. with Fig. [2]) and that the main effect of a larger value of £ is that of 
pushing the onset of the vacuum-dominated expansion to larger values of r, with r accc i oc a/£ 
(compare Fig. [6] with Fig. [3]). It is also instructive to contrast the behavior of the vacuum 
variable x(t) in Fig. H] and that of Fig. |5l where the latter figure displays the "time-lag 
effect" because of the finite dissipative coupling constant C (cf. the heuristic discussion of 
the paragraph starting a few lines under Eq. (4.5) in Ref. |6[). 

As the model evolution for r > Tf reez0 = 3 is perfectly standard (described by an FRW 
universe with ponderable matter and an effective cosmological constant), the crucial segment 
of the numerical solution is over the interval [r min , Tf reeZ e]- The numerical data for the function 
values are given in Table Ull where the £ = oo values refer to the solution of the equations 
(I3.16P derived in Sec. IIII Dl The functions from Table UTI are observed to converge for £ > 10 4 . 
The convergence of the vacuum energy density results is also shown by Fig. [71 With ry (rf reeZ e) 
values for parameters £ both below and above the "realistic" number £ = 10 60 , the following 
estimate is obtained by interpolation: 

£=10 60 

rv(r frceze ) « 0.051, (4.2) 

Table [II] 
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for the model parameters and boundary conditions mentioned in the caption of Table [TT1 

The function values at r freeze calculated from the £ -1 = ODEs (13 . 16[) are given in 
Table HTT1 for a wide range of values of the initial ratio [rMi/ r M2\ at r = r min = 0.1. These 
function values can be expected to approximate the physical (£ _1 ~ 10~ 60 ) values with an 
accuracy of one per mill or better, at least, for initial ratios [rMi/rMi] of order unity. 

As far as the dimensionless remnant vacuum energy density r v (r {reeze ) is concerned, the 
values quoted in Table II III constitute the complete solution of the problem where the initial 
ratio [r mi I f M2\{.Tmm) controls the relative strength of the kick at t ~ t ew , assuming the 
validity of the phenomenological ODEs (13. 7p for coupling constant ( = 0(1) and taking 
model functions k m1 (t) and j(t) from (13.131) and (13.141) . respectively. The dimensionful 
remnant vacuum energy density py (tf Teeze ) requires knowledge of the absolute energy scale 
E ew used in the rescaling of the variables, as will be discussed further in the next section. 



V. DISCUSSION 
A. General case 

In the scenario considered [H,^, the theoretical value of the effective cosmological constant 
(remnant vacuum energy density) is given by 

A thcor ^ = lim $"*(t) = r- m (£ cw ) 8 /(£ Planck ) 4 , (5.1) 

t— >oo 

with ry um a number obtained by numerically solving ODEs of the type of (13.71) . Equating 
the theoretical value (15. ip with the experimental value A cxp w (2 meV) 4 from observational 
cosmology the following estimate for the required energy scale is obtained: 

E m = (A-/ rr ») V ( E^y* „ 3.2 TeV ( p ) "* ( ) "* . (5.2) 

For the moment, the numerical value r™ m = 0.051 in (15. 2p is purely for illustrative purpose. 

Clearly, the calculation of the present paper relies on many assumptions, but it appears 
that values of the order of unity for the dimensionless energy density ry um are quite reason- 
able. The value r y um = 1.0 corresponds to E cw « 2.2 TeV, according to (15. 2p . On the other 
hand, values r™ 111 ^> 1 appear unlikely, at least, in the present framework [see the discussion 
in Sec. |IIIE| leading up to (13.201) ]. Note that an Ty um value of order 10 10 would be required 
in (15. 2p to bring E ew down to the order of magnitude of standard-model particle masses, 
m SM = 10 2 GeV. 

Taking r y um < 1 for granted, the correct reading of (I5.2p is that of a lower bound, 

£ ew >2TeV, (5.3) 

since, without further input, the value of r™ m can be made arbitrarily small (for example, 
by taking a sufficiently small value for the initial energy density tmx m Table HTT1) . Indeed, 
the main uncertainty (apart from the unknown E ew value) is the dynamic importance at 
t ~ t cw of the nonrelativistic (M ~ E ew ) type-1 particles compared to that of the relativistic 
type-2 particles. 
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B. Special case 



In order to get further predictions, the following three assumptions can be made. First, 
assume the type-1 and type-2 particles to have been ultrarelativistic and in thermal equi- 
librium for T ^> E cw , so that their energy densities are given by 

P Mn = (?r/30)N eS , n T\ (5.4a) 

with bosons (b) and fermions (/) of particle type n — 1, 2 contributing as follows: 

N ea , n = 9n,b + (7/8) 9n,f , (5.4b) 

b f 

in terms of the numbers of degrees of freedom of the particles (gj, = 2 for the photon). Then, 
the relevant energy- density ratio Pmi/pm2 before the kick starts is simply given by the ratio 
of the respective effective numbers of degrees of freedom, N e ^i/N e s >2 . 

Second, assume the type-1 particles to have approximately the same mass and a mass 
scale 

M ~ E ew . (5.4c) 

As discussed in Sec. IIII A[ these type-1 particles can be thought to consist of a mix of 
different unstable particles. What matters here, though, is their average thermodynamic 
properties as given by the prescribed EOS function kmi from A3. 13[) . See App. [A] for a 
realistic setup. 

Third, assume the massless type-2 particles to correspond to those of the standard model 
(m < m S M < E ew ), so that 

iV eff , 2 = N cS , SM = 427/4 ~ 10 2 . (5.4d) 



See, e.g., Ref. [20| for the count of the degrees of freedom in the standard model. 

With these assumptions, there are only two unknowns: the numerical value of the energy 
scale E ew and the effective number iV e fr,i of type-1 degrees of freedom. The first quantity, 
E cw , sets, according to A3. 1 j) . the physical time t ~ t cw when the 'kick' of the vacuum energy 
density occurs (the kick mechanism p relies on the change of the Hubble expansion rate by 
type-1 mass effects). The second quantity, N e s,i, controls the initial energy- density ratio: 

[rA4i/r M 2}(r min ) = N cSA /N eS , 2 ~ N cSA /10 2 , (5.5) 

where the dimensionless cosmic time r m i n is taken before the kick starts and r denotes the 
dimensionless energy density according to (I3.12p . 

For a substantial number iV eff x = 10 2 of these ultramassive type-1 particles (possibly 



corresponding to partners of the standard-model particles from broken supersymmetry 21 



22 



), the initial energy-density ratio is given by [rMi/?"M2] ("Train) ~ 1 an d the dimensionless 
remnant vacuum energy density ry um is found to be of order 5.1 x 10~ 2 [see Table IIII] . 
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This particular ry um value requires, according to (15. 5J) , an £ cw value of order 3.2 TeV, in 
order to reproduce the experimental value of the cosmological constant. A similar E ew 
value is obtained in App. [A] if the prescribed (artificial) EOS function (13.131) is replaced 
by a physically-motivated EOS function. Only for iVeff,i = 1 (corresponding to a single 
ultramassive real scalar) does the remnant vacuum energy density drop to such a low value, 
r™ m ~ 2.0 x 10 -5 [Table ITO] . that the required energy scale becomes significantly larger, 
E ew -8.5 TeV. 

Hence, if all particles have initially been in thermal equilibrium and if the type-1 particles 
with mass scale M ~ E ew have an effective number of degrees of freedom iV eff)1 > 1, the 
required energy scale E ew from (15. 2p lies in the following range: 

„ I prescribed kick , , 

Eew k ff .i>i, ^ cff . 2 =io 2 ~ 3 " 9 TeV ' ( 5 - 6 ) 

assuming a dissipative coupling constant ( of order unity. The trend is, as expected from 
( 15. ip . that a smaller number iV effjl requires a larger energy scale E ew . Table ITVl presents 
the required energy scales E ew for selected values of N e g t i, with the understanding that the 
quoted numbers for E ew are only indicative because of the many assumptions made along the 
way (some of which may be more reasonable than others). See App. [A] for further discussion 
of some of the systematic uncertainties involved (its very last table complements Table [TV] 
of this section). 



C. Outlook 

If the observed "cosmological constant" results from dynamics at cosmic temperatures of 
order M ~ E ew , then some set of differential equations must be relevant. In the framework 
of g-theory, a particular set of differential equations has been proposed. It appears to be 
impossible to have very much simpler differential equations which achieve the same result 
and the phenomenological equations used here can be expected to carry some of the essential 
ingredients. If so, the estimates from Table HVI suggest the need for new physics with particle 
masses at the TeV-scale. 

Particle-collider experiments are called upon to confirm or exclude the existence of these 
TeV-scale particles and, if confirmed, to determine their characteristics. Knowing the char- 
acteristics of the new TeV-scale particles (assuming their detection), the main task for 
theorists would be to derive the relevant particle-production effects contained in the simple 
phenomenological equations considered here or to find the appropriate generalizations of 
these equations. 
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Appendix A: Dynamic kick 
1. Introduction 

The goal of this appendix is to present a calculation for the generation of the remnant vac- 
uum energy density by a more or less realistic kick from dynamically generated ultramassive 
and unstable type-1 particles. The description of this dynamic kick is rather involved, but 
the final ODEs are similar to those of the main text for a prescribed kick, with the crucial 
role again being played by the phenomenological dissipation function 7 (the other functions 
entering the ODEs will now be given different notations, for clarity). The heuristics of these 
ODEs will be discussed in the last paragraph of Sec. IA4I 



2. Mass spectrum and EOS functions 

The massless type-2 particles from the model introduced in Sec. IIII Al are considered to 
correspond to those of the standard model (mass scale msM <C E ew ) and the round number 
N 2 = iV"eff2 = 10 2 will be used. The ultramassive type-1 particles are considered to arise 



from broken supersymmetry [21|, |22|, so that Ni = iV effl = iV effj 2 = 10 , and they can 



be taken to be bosons (the standard- model particles being mostly fermions (20j). Just as 
discussed in Sec. IV B\ the type-1 and type-2 particles are assumed to have been in thermal 
equilibrium before the generation of the vacuum energy density starts. 

A general type-1 mass spectrum is given by the effective numbers rij > of particles 
with dimensionless masses m; = Mi/E ew , for which the following two constraints hold: 

n i = N u (Ala) 

i=a,b,c,... 

JjT U i m i = l l ( Alb ) 



where the last constraint ensures that the average dimensionful mass equals E ew . 

For simplicity, consider two cases: case A with two different mass values and case B with 
a single mass value, by definition equal to E ew . (The generalization to a general type-1 mass 
spectrum will be obvious.) The specific numbers are chosen as follows: 

case A: (n a , m a \ n&, rrib) = ( 40, 2; 60, 1/3), (A2a) 

case B: (n a , m a ; n b , m h ) = (A2b) 

The case-B partition {n a , rib} of 100 is arbitrary, at least, for the dynamic ODEs considered 
here. It will be seen, later on, that case A and case B give more or less the same rem- 
nant vacuum energy density, which even holds for the extreme version of case-A having 50 
particles with m a = 2 and 50 with nib — 0. 

Next, take, instead of the prescribed EOS function ( 13 . 1 3[) used in the main text, the exact 
EOS function or, at least, a controlled approximation of it. In fact, the following rational 
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function of the variable 9i = T/Mi will be employed for type-1 subspecies i G {a, b, c, . . .}: 



K M ii(0i) = 1-3 w Mli (Pi) , (A3a) 

"-"W = 3 / + + g ** +g . (A3b) 

(a, /?) = (0.625, 1.91) . (A3c) 

Expression ( 1A3bl) with constants ( lA3cj) gives an accurate approximation (better than 1 per 
mill) to the exact equilibrium result for the EOS parameter wmu = Pmu/pmu of bosons 
at zero chemical potential (see, e.g., Sec. 3.3 of Ref. jl9|). For the EOS parameter Wmu of 
fermions, the constants would be (a, (3) = (0.770, 2.15), but the approximation would be 
somewhat less accurate (still better than 5 per mill). 

As the type-2 particles are massless and adiabatically expanding in the early phase (for 
t < TmiTi <C 1), their initial energy density can be used to normalize the temperature-over- 
mass variable 9i used in ( 1A3al) . 

0i = ET" —T-r PM2 tmin 77^ = — T M 2 T min , A3d 

where rrii E ew is the physical mass Mi of the type-1 subspecies considered. As the scale factor 
a(r) evolves with dimensionless cosmic time r, so do d{ and kmu in the above equations. 
Recall that the scale factor a(r) is defined by h = a/ a and that the temperature of a 



noninteracting gas of "photons" drops as T(r) oc 1/ a(r) in the expanding FRW universe 19 



3. Model functions 

Because the EOS function ~kmu does not vanish for r > Tf reeze , the appropriately modified 
Maxwell equation needs to be used, which has already been discussed in Ftn. |2j But kmu 
also does not vanish exactly in the early phase (as long as a > or T < oo). In order to 
have a standard radiation-dominated FRW universe in the early phase before i ew , the initial 
density of type-1 particles must be taken to be strictly zero. The further assumption is, then, 
that type-1 particles are generated dynamically just before t ew and thermalized rapidly, so 
that the type-li energy density relative to that of the type-2 particles is approximately 
given by rii/N 2 before Kmu starts to differ significantly from 0. 

For the ODEs to be given shortly, the rapid energy transfer from type-2 particles to 
type-1 particles is governed by a coupling constant A21 > and a burst function cD(r), 
which is taken to peak around r = T21 1. For simplicity, define this burst function Q 
using the previous function 7 from (13.14bl) for a fixed value of Tf reeze : 

w(r) = T2i(t)[1-72i(t)], (A4a) 
7 21 (r) = 7(r/r 2 i), (A4b) 
lij) ee 7 (r) I (A4c) 

1 'freeze 
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The subsequent decay of type-1 particles into type-2 particles is governed by a coupling 
constant A12 > and a decay function v{r). For simplicity, this function is taken to have 
the following form: 

u( T ) = 1 - 7(r) . (A4d) 

The hats on the above functions indicate that they are independent of the value of the dis- 
sipative coupling constant (, whereas 7 from (13.141) does depend on it. The model functions 
used will also be shown in the plots of the numerical results later on. 

4. ODEs and boundary conditions 

With a type-1 mass spectrum of the form (IA2j) and the dynamic EOS function Kmh from 



(1A3I) . the dimensionless ODEs are now taken to be the following: 

3 (h + 2h 2 ) = 7 x + (1 - 7) (V 1 x 2 + \ K M ia r M \a + \ Kmu r M ib) , (A5a) 

TMla + (4 - «Aflo) h TMXa = (Nia/Nx) (\ 2 l U>r M 2 ~ (C/t) [ x ]) ~ ^12 Vr M la , (A5b) 

rMib + (4 - Kmu) hr Mlb = (Nib/Nx) ^A 2 i cur M2 - ((/y) [x]j - \iivr M ib , (A5c) 

r M2 + 4:hr M 2 = -A 2 i ur M 2 + Ai 2 v(r M ia + TMib) , (A5d) 

3 (r 1 hx + h 2 ) = r 1 x 2 /2 + r Mla + r Mlb + r M2 , (A5e) 



where the dissipation function 7 is given by (13. 14ft and the other model functions by (1A4I) . 
The basic structure of these ODEs is identical to that of (I3.15P with a prescribed EOS 
function A3. 13|) . 

As the physical hierarchy parameter £ can be expected to be very large (perhaps of order 
10 60 ), the £ _1 = ODEs are especially relevant: 

3 (h + 2h 2 ) = 7X + (1-7)U K Mla [3 h 2 - r M2 - r Mlb ) + \ H Mlh r M \bj , (A6a) 

x = -C 1 h 7 2 (2 x - K M \a [3 h 2 - r M2 - r M u] ~ ^mu r M u ) , (A6b) 

r M ib + (4 - kmu) h r M \b = (N lb /N 1 )(\ 21 ur M2 - - ^i2vr M u, ( A6 c) 

r M 2 + ^hr M 2 = -A 2 iSr M2 + \i 2 v [3h 2 - r M2 ] , (A6d) 

with rj\/ la following from the solution of these ODEs by the corresponding Friedmann equa- 
tion, 

TMXa = 3 h 2 - r M 2 - r M u ■ (A7) 
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For both sets of ODEs, (1A5j) and (1A6j) . the following boundary conditions on h, Tmu, 
Tm2, and x are taken to hold at a time r = r min <C 1: 

h(T min ) = 1/2 (rnon)" 1 , (A8a) 

nVfliOmin) = 0, (A8b) 
»"M2(Tmm) = 3 [/i(r min )] 2 , (A8c) 

x(r min ) = 0, (A8d) 

which correspond to a standard radiation-dominated FRW universe. It must be emphasized 
that there is essentially no free parameter in these boundary conditions, the precise value of 
Tmin being irrelevant as long as it is sufficiently small compared to 1 (so that t min <C £ ew ). 

Before turning to the numerical solutions, it may be helpful to discuss very briefly the 
heuristics of the ODEs (jA5j) . the discussion of Eqs. (IA6[) and ( TA7J) being similar. Starting 



from the values (1A8I) . the A21 oj terms in (lA5bj) and (1A5cl) generate a nonzero value of the 
type-1 matter energy density TMia + ^Mib, which then gives a nonzero value of the vacuum 
variable x from (lA5ep . The values of rMu(r) and x(r) peak around r = 1, which provides 
a dynamic version of the electroweak kick, improving upon the prescribed (artificial) kick 
discussed in Ref. j^] and the main text of the present article. 6 Finally, the value of x 
approaches a constant positive value and so does the vacuum energy density ry = \x 2 1 
whereas Tmu + r A/i& drops to zero because of the decay terms A12 v in ( 1A5bl) and ( 1A5cl) . 



5. Numerical solution 

Numerical results for the case-A mass spectrum (lA2aj) and hierarchy parameter £ = 
10 2 are given in Fig. [8] for the standard (nondissipative) g-theory ODEs with ( = and 
in Figs. [9HT0l for the modified (dissipative) g-theory ODEs with ( = 2. [The coupling 
constant value A 2 i = 18 has been chosen to give r Mla + r MU ~ r M2 at r = 0.25, as would 
approximately correspond to thermal equilibrium. The value of A 12 = 2 has been chosen so 
that most type-1 particles have decayed by a cosmic time r = 3. But the results for ry um 
are more or less constant for A12 < 20, which includes the nonrealistic case A12 = of stable 
(massive) type-1 particles.] Figures M and [TD] are seen to give a somewhat lower value of 
the dimensionless remnant vacuum energy density r™ 111 = ry(r freez e) than Figs. [2] and [3j one 
reason being that the maximum value of kmi t m\ is less for the dynamic kick considered 
than for the prescribed kick (cf. Ftn. [6]). 

Numerical results for the case-A mass spectrum (1A2al) and the physically relevant hi- 
erarchy parameter £ = 00 are given in Fig. [TTJ with ry um of order 1.34 x 10 -2 for ( = 2. 



6 It is clear from (|A6b[) . in particular, that the kick of the vacuum energy density is driven by the combination 

Y^i^Mli r Mli- 
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Similar results (Fig. [T2~]) are obtained for the case-B mass spectrum (lA2bl) . with ry um of 
order 1.10 x 10~ 2 . As mentioned in Sec. fl] Figs. [Til and [T2l present the best (i.e., most 
realistic) calculations of this article. 

According to ( I5.2|) . an average value r™ m ~ 1.2 x 10~ 2 has a required E ew value of 
order 3.8 TeV, which is 20% above the value 3.2 TeV from Table IIVI The mass values 
corresponding to E ew ~ 3.8 TeV are M la ~ 7.7 TeV and M lb ~ 1.3 TeV for the case-A 
spectrum and M\ a ~ 3.8 TeV for the case-B spectrum. With the dynamic kick as defined in 
this appendix, the obtained value E ew ~ 3.8 TeV depends primarily on the assumption that 
the dissipation function 7 enters the ODEs ( 1A6|) in the way shown and that the relevant 
coupling constant ( is of order unity. 

Consider, then, the ( dependence of the predicted E ew value. As explained in the fourth 
paragraph of Sec. IIII E[ the value of the remnant vacuum energy density r™ m (() can be 
expected to drop to zero for C i and £ — > 00. From the numerical solution of ODEs ( IA6P 
for the case-A mass spectrum flA2al) and the same boundary conditions as in Fig. [Til the 
following rough fit of r^ m (Q for £ > 0.02 is found: ry(Q = £ 5 / (a( 7 + /?) with constants 
(5, 3D = (13.29, 680.3). More precisely, the r™ m (0 val ue peaks at 1.34 x 10~ 2 for ( = 2 
and drops to 2.46 x 10" 7 for ( = 0.2 and to 5.82 x 10~ 4 for C = 20. [Note that, for ( < 0.2, 
the time scale Tf reeze from f )3.14ap may become unreasonably large and that, for ( 3> 20, the 
particle-production effects on the right-hand side of flA6cj) . for given value of x oc pv, may 
be larger than can be expected from weakly-coupled particles.] The corresponding results 
for E ew are given in Table |V] which, together with Table IIV[ gives an idea of the systematic 
uncertainties involved (compare the E ew values for N cS)1 = 10 2 and ( = 2 from both tables). 

For the dynamic kick considered in this Appendix, the estimate is thus 

„ I dynamic kick a -, ^ m i r i a /~,\ 

E ew J v 1n2 a ~ 4 - 15 TeV , ( A9 ) 

I "off,i=-Voff,2=10 i case— A y ' 

assuming that the dissipation function 7 enters the ODEs flA6j) in the way shown and that 
the coupling constant ( lies between 0.2 and 20. The corresponding case-A mass values are 
Mi a = 2 -E ew and = 1/3 E ew . The case-B mass spectrum (lA2bj) . with a unique mass 
value Mi a = -^ew, gives the same E ew estimate as in (IA9[) . 
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TABLE I: Function values of h(r), r^i(r), tm2(t), and ry(r) at r = Tf ree ze = 3 from the numerical 
solution of the ODEs A3. 15f) for dissipative coupling constant £ = 2, hierarchy parameter £ = 10 2 , 
and various values of the boundary condition rjvn/ r M2 at r = r min = 0.1. The other parameters 
and boundary conditions are given in the caption of Fig. [TJ The numerical accuracy of the quoted 
function values is estimated to be equal to ±1 in the last digit. 





TMxItM^ (V m i n ) 


h 




rM2 


r v 


10 2 


io- 1 


0.1702 


0.01329 


0.07357 


0.001473 


10 2 


1 


0.1849 


0.06321 


0.03893 


0.04407 


10 2 


10 1 


0.1983 


0.1007 


0.01577 


0.1441 


10 2 


10 2 


0.2008 


0.1070 


0.01220 


0.1705 



TABLE II: Same as Table HI but now for a fixed boundary condition [r mi / r M2\{tmm) = 1 and 
various hierarchy parameters ^ ranging from 1 to 10 6 . The entry for £ = oo has been calculated 
from the ODEs (|3.16p and the algebraic equation (I3.17j) . 





f"Ml/rM2] (Tmin) 


h 


TMl 


TM2 


rv 


i 




1 




0.1758 


0.04279 


0.04394 


0.006026 


10 1 




1 




0.1817 


0.05582 


0.04069 


0.02584 


10 2 




1 




0.1849 


0.06321 


0.03893 


0.04407 


10 4 




1 




0.1859 


0.06530 


0.03841 


0.05052 


10 6 




1 




0.1859 


0.06533 


0.03841 


0.05061 


oo 




1 




0.1860 


0.06533 


0.03841 


0.05061 


TABLE III: Same as Table ED but now from the £ _1 = ODEs f|3.16|) and Eq. (I3.17p. 


e 


Tmx/tmz] (Tmin) 




h 






ry 


oo 




10~ 3 




0.1667 


0.0001577 


0.08322 


2.058 x IO -7 


oo 




IO" 2 




0.1671 


0.001557 


0.08219 


2.021 x 10~ 5 


oo 




io- 1 




0.1704 


0.01384 


0.07330 


0.001699 


oo 




1 




0.1860 


0.06533 


0.03841 


0.05061 


oo 




10 1 




0.1995 


0.1036 


0.01583 


0.1647 


oo 




10 2 




0.2020 


0.1100 


0.01243 


0.1948 


oo 




10 3 




0.2023 


0.1107 


0.01207 


0.1983 
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TABLE IV: Estimates for the energy scale E ew of type— 1 particles with mass scale M ~ E ew 
as a function of their effective number of degrees of freedom iVeffi- The inverse of the hierarchy 
parameter (|3.1ip is taken to vanish, £ -1 = 0, and the dissipative coupling constant is assumed to 
have a value of order unity, specifically £ = 2. The kick from the type— 1 particles is modeled by the 
prescribed EOS function kmi(t) from f)3. 13j) with parameters listed in the caption of Fig. [TJ Both 
ultramassive type- 1 and massless type-2 particles are assumed to have been in thermal equilibrium 
for T 3> E ew and the effective number of type-2 particles is taken as iV e ff 2 = 10 2 . The -Eew estimates 
are obtained from (|5.2p . using the ry values from Table UTTl for [tmi /tm'A (jmm 

) = jV effil /10 2 and 

taking the measured value of the cosmological constant to be A exp = (2 meV) 4 . 



c 




^ prescribed kiek ^yj 


2 


1 


8.5 


2 


10 1 


4.9 


2 


10 2 


3.2 


2 


10 3 


2.8 


2 


10 4 


2.7 



TABLE V: Same as Table HV] but now for the dynamic EOS function kmi( t ) of Appendix lAl 
The total effective numbers of degrees of freedom are iV^i = Ngg^ = 10 2 and the type-1 mass 
spectrum is given by (| A2aj) . The estimates for the energy scale E ew are given for = and 
three values of the dissipative coupling constant £, the other coupling constants being listed in the 
caption of Fig. [9j 
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N eS ,i 


7-t dynamic kick rrp t n 

E c 4 [TeVJ 
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FIG. 1: Numerical solution of the ODEs (13.150 with a prescribed EOS function kmi, for vanishing 
dissipative coupling constant C = and trivial function 7(7*) = 1. Model parameters are {£, A12} 
= {100, 8/100} and the model function kmi(t) is defined by f)3. 13|) with parameters {k c , t c , At} 
= {1, 1, 18/10}. The ODEs are solved over the interval [r m i n , r max ] = [0.1, 3] with the following 
boundary conditions from (|4.ip at r = T m \ n = 0.1: {x, h, a, rjyi, 7^2} = {0, 5, 0.3, 37.5, 37.5}. 
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FIG. 2: Same as Fig. [H but now for nonvanishing dissipative coupling constant £ = 2 and nontrivial 
function 7(7") defined by (|3.14|) with parameter Too = 2, giving Tf ree ze = 3. 
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FIG. 3: Same as Fig. [2l but also evolved for r > Tf reeze = 3 with the standard FRW equations 
[given by (|3.15|) for x = 0, kmi = 0, and 7 = 0]. 
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FIG. 4: Same as Fig. [H still with ( = but now for a larger hierarchy parameter £ = 10 



0.8 
0.6 
0.4 
0.2 

°C 

2 
1.5 
1 

0.5 

°( 


rh 


0.4 
0.3 
0.2 
0.1 

0.8 
0.6 
0.4 
0.2 


X 


1.2 
0.9 
0.6 
0.3 

1.2 
0.9 
0.6 
0.3 


a 4 tMi 


1 

0.75 
0.50 
0.25 

! T °( 
1 

0.75 
0.50 
0.25 


*M1 


















































































































t 
























i 






X 




















































) 0.5 1 1.5 2 2.5 C 

a 


) 0.5 1 1.5 2 2.5 ; 
10 r v 


) 0.5 1 1.5 2 2.5 C 

a 4 r M2 


) 0.5 1 1.5 2 2.5 3 

7 


































































C 


































































































































) 0.5 1 1.5 2 2.5 : 


) 0.5 1 1.5 2 2.5 C 


) 0.5 1 1.5 2 2.5 C 


) 0.5 1 1.5 2 2.5 3 ' 



FIG. 5: Same as Fig. [21 still with C = 2 but now for a larger hierarchy parameter fj = 10 4 . 
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FIG. 6: Same as Fig. [31 still with £ = 2 but now for a larger hierarchy parameter £ = 10 4 . 
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FIG. 7: Remnant vacuum energy density ry(Tf reeze ) from Table ITT1 plotted against the compact 
parameter r\ = log 10 £/(| log 10 £| +3) G [0, 1) defined in terms of the hierarchy parameter £ G [1, oo) 
of the theory. The realistic value £ = 10 60 from (|3.1ip corresponds to r\ = 20/21 ~ 0.9524. 
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FIG. 8: Numerical solution of the ODEs (|A5p with the dynamic EOS function Kmi from (|A3|) . for 
vanishing dissipative coupling constant £ = and trivial function j(t) = 1. Model parameters are 
{£> ^21) ^12}= {100, 18, 2} and model functions are defined by (|A4p with parameter r 2 i = 1/10 
[the function £5(t) peaks below r = 0.3 and ?(t) rises from to 1]. The case-A type— 1 mass 
spectrum (|A2a|) has been assumed, with type-la functions shown by full curves and type-lb 
functions by dashed curves. The ODEs are solved ov6r th.6 interval [r m i n , T max ] = [0.01, 3] with 
the following boundary conditions from (|A8|) at r = r m \ a = 0.01: {x, h, a, r^/ia, ^Mife, f~M2} = 
{0, 50, 0.1, 0, 0, 7500}. 
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FIG. 9: Same as Fig.El but now for nonvanishing dissipative coupling constant ( = 2 and nontrivial 
function 7(7") defined by (|3.14|) with parameter Too = 2, giving Tf reeze = 3. 
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FIG. 10: Same as Fig. O but also evolved for r > Tf reozo = 3. 
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FIG. 11: Same as Fig. [9l still for the case-A mass spectrum (|A2ap but now solving the £ 1 
ODEs (|A6P and using the algebraic Eq. (|A7p . 
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FIG. 12: Same as Fig. [TT1 still with £ 1 = but now for the case-B mass spectrum (|A2bp . 
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